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Abstract. To simulate the transient enhanced diffusion near the surface or interface, a 
set of equations describing the impurity diffusion and quasichemical reactions of dopant 
atoms and point defects in shallow ion-implanted layers is proposed and analyzed. The 
diffusion equations obtained take into account different charge states of mobile or immo- 
bile species and drift the mobile species in the built-in electric field and field of elastic 
stresses. The absorption of self-interstitials on the surface and drift of the defects due 
to elastic stresses result in the nonuniform distributions of point defects. It was shown 
analytically and by means of numerical calculations that consideration of the nonuniform 
defect distributions enables one to explain the phenomenon of "uphill" impurity diffusion 
near the surface during annealing of ion-implanted layers. The performed calculations of 
the boron concentration profile after annealing of a shallow implanted layer agree well 
q ■ with the experimental data confirming the efficiency of the proposed equations. 
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oo : 1 Introduction 

q 

Using low-energy high-fluence ion implantation with the following rapid thermal anneal- 
ing, one can produce the active regions of silicon devices characterized by very shallow 
junctions (~ 0.1 fiva) and high dopant concentrations [U [2j [3j. During annealing the 
transient enhanced diffusion (TED) occurs. Therefore, the final dopant distribution is 
determined both by ion implantation and transient enhanced diffusion parameters. As 
the lateral dimensions of integrated circuits are scaled down to the submicrometer range, 
the need for accurate modeling of silicon doping is increased. The well-known and widely 
used models of TED during rapid thermal annealing of semiconductor substrates very 
often lead to the results which disagree with the experimental data for low-energy high- 
dose ion implantation, especially in the near surface region. For example, the calculations 
presented in Refs. jU, [2] provide some evidence of disagreement with the experimental 
data in the vicinity of the interface. The difference is mainly due to the use of inade- 
quate clustering models for the description of high concentration diffusion pE], due to the 
influence of interfaces on the defect distributions [5], due to dopant atom trapping by 
immobile sinks [6] , and due to the influence of stresses on the dopant and defect diffusion 
[8] . The principal goal of this investigation is to obtain a set of equations describing 
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impurity diffusion and quasichemical reactions of dopant atoms and point defects near 
the surface or interface to simulate adequately the transient enhanced diffusion in shallow 
implanted layers. 



2 System of equations 

It is proposed that during the initial stage of annealing the ion-implanted silicon sub- 
strates a thin damaged layer is formed near the surface or in the vicinity of the interface. 
For example, such a damaged layer can be formed due to incomplete solid phase recrys- 
tallization of silicon regions heavily doped by As [9] or due to the near surface absorption 
of mobile point defects created by boron implantation. During annealing this damaged 
layer can absorb mobile point defects and impurity atoms, that results in the formation of 
nonuniform defect distributions and electrical deactivation of impurity. Moreover, if the 
damaged layer contains a great amount of defects and impurity atoms, significant stresses 
can arise in the region adjoined to the surface. To describe diffusion of point defects and 
impurity atoms in this intricate system, the following set of equations is proposed: 
1. Expression for the total concentration of impurity atoms C T 



iAD 



c 1 = c + c + c 

where C is the concentration of substitutional^ dissolved impurity atoms, C AC and C AD 
are the concentrations of impurity atoms incorporated into clusters and bound to the 
extended defects, respectively; 

2. Conservation law for impurity atoms incorporated into clusters (precipitates) 
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3. Conservation law for impurity atoms bound to the extended defects 
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4. Equation of diffusion for impurity atoms due to the formation, migration, and 
dissociation of the pairs "impurity atom — point defect" which are in equilibrium with 
the substitutional^ dissolved impurity atoms and point defects [TO]. For the two-stream 
diffusion govern by the vacancies and self-interstitials this equation has the form 
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5. Equation of vacancy diffusion [TT] 
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6. Equation for diffusion of self-interstitials [TT] 
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where the effective coefficients of the equations are presented as follows: 
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Here 5 ,j4c and G AC are the rates of impurity atom absorption due to the cluster forma- 
tion and generation of separate impurity atoms during cluster dissolution, respectively; 
S AD and G AD are the rates of impurity atom absorption due to the formation of extended 
defects and generation of separate impurity atoms during extended defect annealing, re- 
spectively; C Vx and C lx are the concentrations of vacancies and self-interstitials in the 
neutral charge state normalized to the equilibrium concentrations C^ x and C\ x respec- 



tively; D E (x) is the effective diffusivity of impurity atoms due to the vacancy — impurity 
pairs mechanism; D F (x) is the effective diffusivity of impurity atoms due to migration 
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of the pairs "impurity atom — self-interstitials" Df and Df are the intrinsic diffusivi- 
ties; superscripts E and F denote the vacancy and interstitial pair diffusion mechanisms, 
respectively; /3 1 ' F and /3 2 ' are the parameters describing respectively the relative con- 
tributions of singly and doubly charged point defects to the impurity transport; functions 
a E ' F describe the overall influence of heavy doping on the local equilibrium between dif- 
fusing pairs, point defects and substitutionally dissolved dopant atoms and also between 
the point defects in different charge states; v% and vf are the projections onto the x -axis 
of the effective drift velocities of dopant atoms in the field of elastic stresses; \ is the con- 
centration of charge carriers normalized to the intrinsic carrier concentration rif, C B is the 
concentration of impurity with opposite type of conductivity; z cl and m are the cluster 
charge in units of the elementary charge and the number of doping atoms incorporated 
in the cluster; dX' 1 (x) an d are respectively the effective and intrinsic diffusivities of 
point defects; superscripts V and / denote the vacancies and self-interstitials, respectively; 

and t^' 1 are the effective recombination coefficient and the average lifetime of the 
corresponding point defects, respectively; functions a > describe the influence of heavy 
doping on the local equilibrium between point defects in different charge states; v£ and 
v x are the projections onto the x -axis of the effective drift velocities of vacancies and 
self-interstitials in the field of elastic stresses, respectively; £{ and l\ are the average mi- 
gration lengths of vacancies and self-interstitials, respectively; gf and g\ are respectively 
the generation rates of vacancies and self-interstitials. 

As follows from Ref . jl] , to describe the impurity clustering during transient enhanced 
diffusion of As and P, one can use the expression 

C AC = KC DX m ~ zCl C m (14) 

instead of Eq.(j2J), if the annealing duration is not very short (for example, annealing is 
longer than 1 s at a temperature of 1000 °C and higher). A similar expression can be 
used for the description of phosphorus clustering [12J. 

Here K is the constant of the cluster formation reaction; m is the number of impurity 
atoms incorporated in cluster; Cr> is the relative concentration of nonequilibrium defects 
participating in clustering 

(C*Z)l) mi ~ (jDlx _ qD2x 

6d = JS^ • dm = c^ • dm = W • (15) 

where C Dlx and C D2x are the concentrations of point defects and D£ in the neutral 
charge states facilitating the cluster formation and generated during clustering, respec- 
tively; C^ lx and C^ 2x are the equilibrium concentrations of these defects; mi and m 2 
are respectively the numbers of point defects Di and D 2 participating in the cluster for- 
mation. Due to multiplier Cp in the expression (Tl4"|) . it is possible to take into account 
the influence exerted by the nonuniform distributions of nonequilibrium point defects on 
the clustering process and influence of the defect generation during clustering on tran- 
sient enhanced diffusion. Thus, the expression ( Tl4l enables one to consider the coupled 
phenomenon of clustering and transient enhanced diffusion. 
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3 Analysis of equations 



The set of Eqs.flT]) — (Q can describe different processes of diffusion and chemical reactions 
in semiconductors including "uphill" impurity diffusion in the vicinity of the interface. The 
proposed equations take into account different charge states of all mobile and immobile 
species and also the drift of the mobile species in the built-in electric field and in the field of 
elastic stresses, although only the concentrations of neutral defects are given in the explicit 
form in these equations. The phenomenon of "uphill" diffusion can arise due to the direct 
absorption of impurity atoms by extended defects [6] or due to clustering in the near 
surface region, if these processes are nonequilibrium and proceed in the forward direction. 
To describe this phenomenon one can use Eqs.(j2]) and ([3]), respectively. Moreover, as 
can be seen from Eq.(Tj0), the "uphill" impurity diffusion occurs due to the influence of 
stresses on the pair migration [10J or due to the formation of nonuniform distributions 
of point defects in the neutral charge state [13J. As it follows from Eqs.Q and (jBJ), 
the nonuniform distributions of vacancies and self-interstitials can be formed due to the 
absorption of point defects by the surface (interface) or by extended defects. Moreover, 
the nonuniform defects distributions can be formed due to the stress-mediated migration 
of vacancies and self-interstitials. Thus, a multifactor character of the "uphill" impurity 
diffusion in the near surface region follows from analysis of Eqs.fll]) - - ([61), and it is 
difficult to determine which factor is dominant in the formation of the near surface peak 
of the impurity concentration. However, there is an experimental evidence that the silicon 
surface absorbs self-interstitials [5J, which govern the transient enhanced diffusion. It 
means that in any case the formation of nonuniform distributions of point defects due to 
their absorption on the surface provides a certain contribution to the "uphill" impurity 
diffusion. 

Using an analytical approach, let us evaluate the influence of the nonuniform neutral 
defect distribution on the "uphill" diffusion. Consider the case of TED governed by one 
species of nonequilibrium point defects, namely nonequilibrium self-interstitials. Let us 
suppose for simplicity that the concentration of an impurity providing the opposite type 
of conductivity is approximating zero C B ~ and a F ~ 1. Also, suppose that during 
annealing the formation of neutral clusters or precipitates occurs. Then, in the region of 
high impurity concentration Eq. (pE]) takes the form 

The first term in the right-hand side parentheses describes the impurity flux according 
to the Fick's law, whereas the second one describes an additional flux arising due to 
the nonuniform defect distribution. If this flux is opposite to that due to the impurity 
concentration gradient, the phenomenon of "uphill" diffusion is observed. 

For simplicity, we consider the case of long-time annealing. Then, a steady-state 

dC T 

solution of Eq. (16) for t — > oo and ^ ~ can be derived 

2D F ( X ) C 1 ^ +D F ( X ) C^- = A = const . (17) 
Let us consider the structure "silicon-on-insulator" and suppose reflecting boundary 
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conditions for impurity atoms on the surface (x = 0) and at the interface (x = Xr). Then 
A = 0, and Eq. (17) can be written as 



^- = -b(x)C 
ax 



where 
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We can solve Eq.([T 
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using the method of variable separation 
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C (x) = Cr exp 
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where Cr = C (xr) is the impurity concentration at the interface. 

Let us suppose for simplicity that the distribution of self-interstitials is approximately 
described by the solution of Eq.flE]) with Dirichlet boundary conditions and constant equa- 
tion coefficients taking the influence of stresses and generation of point defects in the near 
surface region as negligible 



C lx (x) = sinh" 1 ( ^ 
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where Cg X and (7^ x are the normalized concentrations of point defects at the surface and 
interface, respectively. 

Substituting expression ffl9l into ff2~Tl) . one can obtain 
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Integral calculation in expression (T23]) gives 



C (x) = Cr exp 
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It can be seen from expression (|24j) that the steady-state distribution of impurity 
atoms varies inversely as the square root from the normalized self-interstitial distribution 
C lx (x). It means that for the defect distribution, similar to the distribution of (22) with 
Cg X << (7jj X , i.e. in case of the absorption of self-interstitials at the surface, the "uphill" 
impurity diffusion and peak of impurity concentration adjacent to the surface must be 
observed. It follows from Eqs. ([1]) and (114"]) that this peak will be more pronounced in 
the case of significant clustering of impurity atoms. 
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4 Simulation results 



Results of the previous analytical investigation are confirmed by the numerical simulation 
presented below. As an example, the redistribution of boron atoms implanted into silicon 
at low energy and high fluence is investigated. For comparison, the experimental data 
obtained by Lerch et al. [H] are used. In Ref. [H] boron ions were implanted into n-type 
silicon substrates of the resistivity 10-20 Qm and orientation (100) at the energy of 1 keV 
and fluence Q = 1.0 xlO 15 cm~ 2 . The annealing process was carried out at a temperature 
of 1050 °C with duration of 10 s. The concentration profile of the boron atom distribution 
was measured by SIMS. This concentration profile is presented in Fig. [1] together with 
the profile calculated by means of the set of equations obtained. 

The finite-difference method [15] is used to find a numerical solution of this system 
of equations. The present paper chooses a more stable implicit method to avoid numer- 
ical instability. Following [15], the right-hand side of Eq.(T4|) and two first terms on the 
left-hand side of Eq. ([6]) which describe the diffusion of dopant atoms and self-interstitials, 
respectively, are approximated by a symmetric difference operator of second order accu- 
racy on the space variable x. To solve the system of the obtained nonlinear algebraic 
equations we use the simplest iterative technique substituting the values of the dopant 
concentration and concentration of neutral self-interstitials determined on the previous 
iteration into the coefficients of nonlinear algebraic equations. 

During simulation it was supposed that transient enhanced diffusion of the boron 
implanted into silicon substrate occurred due to the kick-out diffusion mechanism that 
is mathematically equivalent to the pair diffusion mechanism [16]. We have used the 
following parameters to describe ion implantation and diffusion of the dopant atoms and 
point defects: Q = 1.0 xlO 15 cm" 2 ; R p = 0.00512 /j; AR P = 0.0038 fx; Sk = 0.3; n { = 
1.166xl0 7 ^ 3 ; AALax = 1.3xl0- 05 ^ 2 /s; (3{ = 0.7; # = 0; Cf max = C' a * /C 1 ^ =0.04 
a.u.; l\ = 0.02 \i. Here R v is the average projected range of dopant ions; AR P is the 
straggling of the projected range; Sk is the skewness of boron distribution as implanted; 
Di is the intrinsic diffusion coefficient of boron; C^ x = C^ax/ Ceq * s a maximum value 
of the normalized self-interstitial concentration; C^ x is a maximum value of the neutral 
self-interstitial concentration in the region of impurity diffusion; C| x is the concentration 
of neutral self-interstitials at the surface. The value of equilibrium boron diffusivity at 
1050 °C known from literature = 4.098 xlO -6 ^ 2 /s |17j . This means that the time 
average enhancement of boron diffusion due to the nonequilibrium self-interstitials is 
approximately 3.2 fold. 

As can be seen from Fig. [H the measured boron profile is characterized by the for- 
mation of a local peak of the impurity concentration near the surface due to the distinct 
"uphill" diffusion. It is interesting to note that the impurity profiles calculated by means 
of the computer code SSUPREM-IV included in the desktop program (framework) Athena 
which are presented in [14J do not describe this phenomenon. In our opinion, this dis- 
crepancy with the experiment in the near surface region is due to neglecting the interface 
effect on the defect subsystem that is absolutely inappropriate for the models intended 
to simulate shallow doped layers. It was shown in [5] that the semiconductor surface 
absorbs nonequilibrium self-interstitials governing the transient enhanced diffusion of the 
ion-implanted impurity. As can be seen from expression ff24^) and numerical calculations 
presented in Fig. [TJ taking into account the phenomenon of self-interstitial absorption 
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Figure 1: Calculated concentration profile of ion-implanted boron after thermal annealing 
during 10 s at a temperature of 1050 °C. Solid line — calculated distribution of the total 
boron concentration after annealing; dash-dotted line — distribution of boron after ion 
implantation; dotted line — time average distribution of self-interstitial concentration. 
Circles represent the experimental data taken from Lerch et al. [14J. 

enables one to explain the formation of a near surface peak of the impurity concentration 
and provides a good agreement of the calculated impurity profile with the experimental 
data. It is interesting to note that absorption of boron atoms by extended defects in the 
near surface region also can contribute to the "uphill" diffusion. In this case the nonuni- 
form self-interstitial distribution with decreased concentration gradient can be used to 
explain the experimental data. 
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5 Conclusion 



For adequate simulation of the transient enhanced diffusion near the surface or at the 
interface of ion-implanted layers the set of equations describing impurity diffusion and 
quasichemical reactions of the dopant atoms and defects was proposed and analyzed. 
The proposed diffusion equations take into account different charge states of mobile and 
immobile species and drift of the mobile species in the built-in electric field and the 
field of elastic stresses, although in these equations only the concentrations of neutral 
defects are given in the explicit form. The absorption of self-interstitials on the surface 
and drift of mobile defects due to elastic stresses result in the nonuniform distributions 
of point defects. It was shown, analytically and by means of numerical calculations, 
that consideration of the nonuniform defect distributions makes it possible to explain 
the phenomenon of "uphill" impurity diffusion near the surface during annealing of the 
ion-implanted layers. The performed calculations of the boron concentration profile after 
annealing of a shallow implanted layer agree well with the experimental data including the 
near surface region, where the formation of a concentration peak occurs. The agreement 
with experiment confirms the efficiency of the proposed equations for simulation of shallow 
layer annealing. 
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